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THE STABILITY OP ECCENTRICALLY 
STIFFENED CIRCULAR CYLINDERS 


VOLUME I 
GENERAL 


By 

G. W. Smith and E. E. Spier 
General Dynamics Convair Division 
San Diego, California 


ABSTRACT 

This is the first of six volumes, each bearing the same report number, 
but dealing with separate problem areas concerning the stability of 
eccentrically stiffened circular cylinders. The complete set of documents 
was prepared under NASA Contract NAS8-11181 and furnishes workable design 
and analysis methods for the prediction of instability in such structures. 
The overall multiple-volume report includes design curves, procedures, and 
digital computer programs for the most important buckling modes. The 
material presented in this first volume provides an introductory background 
which puts the subsequent volumes in a proper perspective. 


iii 



GDC- DDG-6 7-006 


TABLE OF CONTENTS 


Section 

Title 

Page 


ACKNOWLEDGEMENTS 

ii 


ABSTRACT 

iii 


LIST OF ILLUSTRATIONS 

vii 


DEFINITION OF SYMBOLS 

ix 


GLOSSARY 

xi 

1 

INTRODUCTION 

1-1 

2 

BUCKLING CRITERIA 

2-1 

3 

BASIC ORTHOTROPIC CYLINDER EQUATION 

3-1 

4 

OVERALL ANALYSIS PROCEDURE 

4-1 

5 

REFERENCES 

5-1 


Preceding page blank 


GENERAL DYNAMICS CONVAIR DIVISION 



GDC- DDG-6 7-006 


Figure 

1 

2 

3 

4 


5 


LIST OF ILLUSTRATIONS 


Title Page 

Panel Instability xii 
General Instability xii 
Anticlastic Bending of a Beam xiv 


Equilibrium Paths For a Perfect 
Isotropic Circular Cylinder Subjected to 
Axial Compression 2-2 

Coordinate System 3-2 


Preceding page blank 


vii 

GENERAL DYNAMICS CONVAIR DIVISION 



GDC-DDG-6 7-006 


DEFINITION OF SYMBOLS 


Symbol 


Definition 


A 11 ,A 22 ,A 12’ A 21 ,A 33 


Elastic constants [see equations 
(3-4)]. 


C 11 ,C 22’ C 12 ,C 21 


D 11 ,D 22 ,D 12 ,D 21 ,D 33 


Eccentricity coupling constants 
[see equations (3-4) J. 

Elastic constants [see equations 
(3-4)]. 


E 

G 

h 

L 


M 

y 


M 


xy 


m 


N, 


THIEL 


(^thiel) 


N 


N 


N 


xy 


Young's modulus. 

Modulus of elasticity in shear* 

Distance between middle surfaces of 
sandwich facings. 

Overall length of cylinder. 

Running bending moment with respect to 
middle surface of basic cylindrical skin 
(acting on sections obtained by passing 
planes normal to the axis of revolution). 

Running bending moment with respect to 
middle surface of basic cylindrical skin 
(acting on sections which lie in radial 
planes) . 

Running twisting moment with respect to 
middle surface of basic cylindrical skin. 

Number of axial half-waves in buckle 
pattern. 

Loading parameter defined in equations 
(3-3), [positive for tensile loading], 

= - , [positive for compressive loading]. 

Applied longitudinal tensile running load 
acting at the centroid of the effective skin- 
stringer combination. 

Circumferential tensile running load acting 
at the centroid of the effective skin-ring 
combination. 

Running shear load acting in the middle 
surface of the basic cylindrical skin. 
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Symbol 

n 

R 


t 

t f 

u,v,w 


x,y,z 

z 


a 



y y 

X 2 * 'yz 


c 

X 


c 

y 

T] 

P 

T1 

s 

V 


DEFINITION OF SYMBOLS 
(Continued) 

Definition 


Number of circumferential full waves 
in buckle pattern. 

Radius to middle surface of basic 
cylindrical skin. 

Thickness of basic cylindrical skin. 

Sandwich facing thickness. 

Reference- surface displacements (see 
Figure 5). 

Coordinates (see Figure 5). 

Parameter defined by equation (3-2). 

Parameter defined in equations (3-3). 

Parameter defined in equations (3-3). 

Parameter defined in equations (3-3). 

In-surface shear strain for middle surface 
of basic cylindrical skin. 

Shear strains in planes normal to the middle 
surface of basic cylindrical skin. 

Longitudinal extensional strain for middle 
surface of basic cylindrical skin (positive 
for pure tension). 

Circumferential extensional strain for middle 
surface of basic cylindrical skin (positive 
for pure tension). 

Parameter defined in equations (3-3). 
Parameter defined in equations (3-3). 
Poisson's ratio. 
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GLOSSARY 

Note: This glossary is meant to apply to all six volumes 

of the report. Separate glossaries are not included 
in each of the individual ' volumes. 

Buckling of Isotropic Skin Panel - The initial buckling of the basic 
cylindrical skin whose boundaries are formed by the longitudinal and/or 
circumferential stiffeners. Buckling of the wall of unstiffened cylinders 
is a special instance of this mode of instability. 

Local Buckling of Longitudinal Stiffener (Stringer) - The initial buckling 
of any leg or arc length of the cross-sectional shape of a longitudinal 
stiffener (stringer). Initial buckling of the outstanding flange of a 
Z-section stringer is an example of this mode of instability. 

Crippling of Longitudinal Stiffener (Stringer) - The final ultimate com- 
pressive failure of a longitudinal stiffener which has a shaped cross section 
and is given sufficient support to prevent panel instability (see definition 
below). The crippling stress is the ultimate average stress for such a 
stringer. 

Panel Instability - This mode of instability manifests itself as a bowing 
of the longitudinal stiffeners (stringers) into one or more axial half- 
waves without any radial displacement of the circumferential stiffeners 
(rings). Hence the axial half-wavelength cannot exceed the spacing between 
rings. Figure 1 depicts the special case where there is one half-wave 
between adjacent rings. Although not frequently encountered in practical 
structures, this mode can involve more than one such half-wave per ring 
spacing. Panel instability may or may not be preceded by buckling of the 
isotropic skin panels and/or local buckling of the stringers. The identifies 
tion “Panel Instability” is somewhat of a misnomer since this terminology 
could easily lead one to the erroneous conclusion that reference is being 
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-Circumferential 
Stiffeners (Rings) 


Axis of Revolution 


Axial 

Load 


Stringer Config- 
uration Just Prior 
to Buckling 


Figure 1 - Panel Instability 

made to the mode which is identified above as "Buckling of Isotropic Skin 
Panel." A more suitable title could be selected but in the interest of 
maintaining consistency with the nomenclature usually found in the literature 
the "Panel Instability" label has been retained in this study* 

General Instability - This mode of instability involves the simultaneous 
radial displacement of both the longitudinal and circumferential stiffeners 
(stringers and rings). As shown in Figure 2, the axial half -wave length 



Axis of Revolution 

Figure 2 - General Instability 
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of the buckle pattern exceeds the spacing between rings. This phenomenon 
may or may not be preceded by buckling of the isotropic skin panels and/or 
local buckling of the stringers. 


Stiffener Eccentricity - The distance from the reference surface to the 
centroid of the appropriate skin-stringer combination. Throughout this 
study, the reference surface has been taken as the middle surface of the 
basic cylindrical skin. Inside versus outside positioning of the stiffeners 
is identified by means of an appropriate sign convention. 


Monocoque - This term comes from the French word meaning "shell only" and 
is used here to identify those configurations which do not incorporate any 
stiffening members (integral or non-integral). Note, however, that a 
monocoque configuration can have orthotropic properties. 

Knock-down Factor - An empirical correction factor which is used in the 
reduction of classical small-deflection predictions to safe design levels. 
This factor is introduced primarily to allow for the detrimental effects 
from initial imperfections. 


Stress Resultants - The six quantities N , N , N , N , Q , and Q 
_r_l " y xy yx x y 

obtained by integration of the infinitesimal loads over the shell wall 

(including skin and stiffeners), and the four quantities M , M , M , and 

x y xy 

M obtained by integration over the shell wall (including skin and 

yx 

stiffeners) of the infinitesimal moments with respect to a selected surface. 
The force stress resultants are expressed in units of force per unit length 
(lbs/in for example) while the moment stress resultants are expressed in 
units of moment per unit length (in- lbs/in for example). 


Shell (or Shell Wall) - Throughout the several volumes of this report, 
repeated use is made of the terms "shell" and "shell wall". These terms 
are used interchangeably and are not meant to refer only to the basic 
cylindrical skin of the stiffened structure. They refer to the effective 
skin-stiffener combination. Whenever it is desired that reference be made 
solely to the basic monocoque cylinder which the stiffeners augment, the 
word "skin" will actually be included in the identification. 
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Anticlastic Bending - Banding into a deflaetad shape for which the 
principal radii of curvature have opposite signs. Bending of an initially 
flat plate into a saddle shape is an eauusple. In addition, for beans, the 
Poisson-ratio effect results in antielastic* bending as depicted in Figure 3 




Figure 3 - Anticlastic Bending of a Beam 
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SECTION 1 
INTRODUCTION 

This is the first of six volumes, each bearing the same report number, 
but dealing with separate problem areas concerning the stability of 
eccentrically stiffened circular cylinders. The areas treated by the 
complete set of documents can be identified from the following listing 
of titles: 


VOLUME 

I - 

GENERAL 

VOLUME 

II - 

BUCKLING OF CURVED ISOTROPIC SKIN PANELS; 
AXIAL COMPRESSION 

VOLUME 

III - 

BUCKLING OF LONGITUDINALLY STIFFENED CYLINDERS; 
AXIAL COMPRESSION 

VOLUME 

IV - 

GENERAL INSTABILITY OF CYLINDERS HAVING 
LONGITUDINAL AND CIRCUMFERENTIAL STIFFENERS; 
AXIAL COMPRESSION 

VOLUME 

V - 

EFFECTS OF INITIAL IMPERFECTIONS; AXIAL 
COMPRESSION AND PURE BENDING 

VOLUME 

VI - 

INTERACTION BEHAVIOR 


The material presented in this first volume is primarily intended to 
provide an introductory background which places the significance of the 
subsequent volumes in a proper perspective. 
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SECTION 2 
BUCKLING CRITERIA 

In the solution of buckling problems, a number of different approaches 
may be taken. The two most commonly used techniques are the minimum 
energy method and the bifurcation concept. The former is based upon the 
theorem of minimum total potential energy which may be stated as follows: 

A conservative system is in a configuration of stable 
equilibrium if, and only if, the value of the total 
potential energy is a relative minimum. 

To apply this theorem, one must formulate the total potential energy of the 
system, impose the mathematical artifice known as a virtual displacement, 
and examine the sign of the second-order energy changes (second variation). 

The second variation must be positive definite (positive regardless of the 
sign of the virtual displacement) for stability to exist. The value of the 
applied load at which the second variation first ceases to be positive 
definite is the critical load for the system. 

The bifurcation concept originally developed by Poincarfe [credited in 
Ref. l] in 1883 constitutes an equilibrium approach to the problem of 
buckling. Any point at which a single equilibrium path branches into two 
or more equilibrium paths is known as a bifurcation point. An example of 
this phenomenon is shown in Figure 4. This figure depicts the equilibrium 
paths for a perfect isotropic circular cylinder subjected to axial compression. 
As a rule, the unbuckled configuration becomes unstable at a bifurcation 
point and the bifurcation approach to buckling analysis simply involves the 
mathematical search for these branching points. In conducting this search, 
one must study the character of the equilibrium behavior. As in the 
solution by Block, Card, and Mikulas [ 2 ], this study may be done with the 


Numbers in brackets [ ] in the text denote references listed in SECTION 5. 
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assistance of energy principles* How ever, such investigations should not 
be misconstrued as constituting a minimum energy approach. In the minimum 
energy method, a so-called second variation is examined. On the other hand. 


Axial 

Compression 

Load 


End Shortening 

Figure 4 - Equilibrium Paths For a Perfect Isotropic 
Circular Cylinder Subjected to Axial 
Compression 

the bifurcation method requires a study only of the so-called first 
variation. That is, in this case the system is subjected to a virtual 
displacement and the first-order change (first variation) in the total 
potential energy is tested for compliance with the principle of station- 
ary potential energy which may be stated as follows: 

A necessary and sufficient condition for the equilibrium 
of an elastic body is that the first-order change in the 
total potential energy of the body be equal to zero for 
any virtual displacement. 
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Although this principle is often used in the application of bifurcation 
theory, it can sometimes lead to a loss of feel for the physical signi- 
ficance of the mathematical operations. Hence, many applications of the 
bifurcation concept are achieved without the use of energy methods. 

Instead, the basic laws of elementary mechanics are used to write the 
governing equilibrium equations for a distorted free-body element. To 
accomplish this, it is only necessary to perform simple summations of 
forces and moments in the several appropriate directions. This, of course, 
is only the first of a series of rather complicated steps in the derivation 
of a final buckling equation but, once having embarked upon this approach, 
each of the subsequent operations likewise retains greater physical trans- 
parency than do the variational techniques. 
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SECTION 3 

BASIC ORTHOTROPIC CYLINDER EQUATION 

The basic orthotropic cylinder equation utilized in Volumes III Csl 
and IV [ 4 ] of this report was developed from bifurcation theory without 
the use of energy concepts. The final equation may be written as follows: 


(^thiel) 


|i + 2n p ^/7 p 2 + y p 4 j 


4a P 


a p 4 (Z) 2 


1 + 2T1 s P 2 + p 4 ] 


(3-1) 


where , 


Z = 


1 - 


c + c 

11 22 


2a tA 22 D 22* ^ P 


12 

2a A 22 (D 22 /Ah) 


1/2 


21 


2a 


(3-2) 


A detailed derivation of this relationship is given in reference 5 where 
the coordinate system shown in figure 5 was used. It might be noted however 
that the final equation developed in reference 5 includes a quantity C , 

Ou 

which is a measure of the shear-center offset from the basic cylindrical 

skin. As written above, the fundamental equation is based on the assumption 

that the stringers and/or rings provide no resistance to y shear de- 

xy 

formations so that C = 0. That is, it has been assumed that all of this 

in-surface type of shear restraint is furnished by the basic cylindrical 

skin. Furthermore, it has also been assumed that transverse shear strains 

Y and Y are everywhere zero. This, of course, is common practice in 
xz yz 

thin-shell theory. However, since this latter assumption is certainly 
not justified for most sandwich-type cylinders, equations (3-1 ) and (3-2) 
must be considered inapplicable to such configurations. 
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Equations (3-l) and (3-2) have been written in their rather compacti 


instructive forms through the incorporation of the following parameters, 
most of which were first proposed by Thielemann [6]: 



T) 


P 


( P 12 + 2 D 3 3 ) 

Vvv 


(3-3) 
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Y 


D 

D 


A 

11 11 
22 A 22 



(3-3) Cont'd. 


The various A. .'s, D. .'s, and C. . *s of equations (3-2) and (3-3) are 
13 13 13 

important fundamental constants which arise out of the following relation- 
ships between the stress resultants and the buckling distortions: 

M x - - D ll 0 - D 12 * c n 5 x + C 12 S y 


a 2 a 2 

M y ■ - D 21 - ”22 ^ - C 21 5 * * c 22» y 

M *y = - 2D 33 tj? 


(3-4) 


C x “ A ll N x + A 12 N y + C 11 


d 2 w a 2 w 

5x 2 + 21 By 2 


A 2 A 2 

, — ,, ow _ ow 

C y = 21 N X + A 22 N y + 12 ^2 + °22 9y 2 


Y = A__ N 
xy 33 xy 

The minus signs in these equations are due to the particular sign conven- 
tions selected in the derivation. In addition it is pointed out that, in 
the above equations, stress resultants with bars above them are centroidal 
values. All other stress resultants are reference surface values. In 
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addition, all of the strain and displacement values are likewise measured 

in the reference surface. The A, . *s and D J , f s of equations (3-2) through 
, . ij ij 

13-4) are usually identified as elastic constants while the C.^’s might be 
suitably described by the terminology "eccentricity coupling constants”. 


It should be noted that equations (3-1) and (3-2) were derived for a 
monocoque orthotropic cylinder. However, these equations will be used in 
subsequent volumes [3,4] for the analysis of discretely stiffened shells. 

The key to success in these applications lies in the means employed to 
evaluate the elastic constants and the eccentricity coupling constants.. From 
equations (3-4) it can be seen that these quantities are dependent upon the 
various structural rigidities of the shell wall. Practical formulas for 
computing these constants are presented in the procedures of Volumes III £3] 
and IV [ 4 ]. However, at this time it is profitable to devote some attention 
to their origins and to examine the formulations which would apply in two very 
special cases. In the first place, it is helpful to note that for an isotropic 
cylinder the elastic constants take on the following forms: 


A - A - — 

11 " 22 ” Et 

A _ A )L 

12 “ A 21 “ Et 


33 


_ 1 _ 

Gt 


D = D 
11 22 


Et‘ 


12 (l-v ) 
l3 


D = D 
12 


vEt 


33 


Gt 

12 


21 

3 


12 (l-v ) 


(3-5) 
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As already noted, the equations of this volume are generally to be considered 

inapplicable to sandwich structures in view of the neglect of transverse 

shear strains (y and y ) in the derivations. However, at this point it 
xz yz r 

is still informative to note that the following formulas could be used to 
find the elastic constants in the very special case of a sandwich configura- 
tion having a core with infinite transverse shear rigidity; 


A 11 

= A = 

22 

2t f E 

A 


V 

12 

" A 2l = ” 

2t f E 

A 33 

= 2^G 


n 

= D 22 = 

Et f h 2 

D n 

2(l-v 2 ) 

n 


vEt f h 2 

°12 

~ °21 “ 

2(l-v 2 ) 

T\ 

Qt f h 2 



D - — - 

33 2 


where 

t f = Facing thickness 

h = Distance between middle- surfaces of 
facings. 

These equations are applicable when the facings are of the same material and 
of equal thickness and this thickness is small compared to h. 

From equations (3-4), (3-5), and (3-6), it should be observed that 

(a) constitutes the reciprocal of the longitudinal extensional 
stiffness per unit length of circumference. 

(b) A ^2 constitutes the reciprocal of the circumferential ex- 
tensional stiffness per unit of axial length. 
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(c) 

(d) 

(e) 

(f) 

(g) 

(h) 


D, 1 constitutes the longitudinal flexural stiffness per 
unit length of circumference. 

D 22 constitutes the circumferential flexural stiffness per 
unit of axial length. 

A 12 and A gl each constitute measures of coupling between 
extensional deformations in the longitudinal and cir- 
cumferential directions. 

D^ 2 and X> 21 each constitute measures of coupling between 
flexural deformations in the longitudinal and circumferential 
directions* 

A 33 constitutes the reciprocal of the in-plane shear stiffness 
of the shell wall. 

Dg 3 constitutes the twisting stiffness of the shell wall* 


In addition it is pointed out that, for applications to skin- stringer-ring 
configurations, the constant is simply the eccentricity (see GLOSSARY) of 
the effective skin-stringer combination. Correspondingly, C 22 is simply the 
eccentricity of the effective skin-ring combination. The C^ 2 and values 
account for Poisson-ratio cross-linking of the eccentricity influences. 

It is helpful to note that equations (3-l) and (3-2) essentially 
comprise an extension to the theory developed in reference 7 for longitudinally 
stiffened circular cylinders. The extension was accomplished in order to 
adapt the solution to circular cylinders having both stringers and rings. 
Equations (3-1 ) and (3-2) can be specialized to the case of a longitudinally 


stiffened circular cylinder by taking C, 


22 


= C =0. 
21 


This gives the result 


/» \ I 1 ^ 2J] ^ *1 * T e 4 ] t _L 

\thiel| c = 4ais o f l + 


a p 4 (Z y 


2T p 2 + 


P 4 


(3-7) 


where 


Z = 


1 - 


11 


12 


2“ U 22 D 22 ,1/!! ^ 


2a A 22 ^ D 22^ A 11^ 


1/2 


(3-8) 
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It is noted that, except for the C 12 term, equations (3-7) and (3-8) are 
identical to equation (36) of reference 7 where the C^ 2 term was evidently 
discarded as a negligible quantity for the' particular test specimens of 
interest there. Further note that when = C 22 = C 12 = C 21 = 0, equations 
(3-l) and (3-2) reduce to the following: 


(^thiel) c = 


[1 + 2Ti p VY P 2 + rP 4 ] 


+ 


<*P 


4 


4aP 


1 + 271(3 + 

s 


e 4 ] 


(3-9) 


This expression is identical to that given by Almroth in references 8 and 9 
except for an apparent typographical error involving the omission of a Vy" 
factor in one term of both references. In addition, it has also been 
established that equation (3-9) is identical to equation (3.2) in a report 
by Appel C 10] . Equation (3-9) can be further simplified by ignoring the 
influences of finite cylinder length, in which case one can arrive at the 
small deflection equation of Thielemann [ 6l which may be expressed as follows 


( s thiel) c = 


1 + 21) VT P 
— . - — p — 


2 + YP 4 


T 1/2 


1 + 21) p 

s 


+ (3- 


(3-10) 


Reference 11 gives a detailed presentation of the mathematical operations 
involved in this simplification. 

To make proper use of equations (3-1) and (3-2), one must realize that 
these relationships simply establish the magnitudes of longitudinal com- 
pressive loads which will maintain an orthotropic cylinder in deflected con- 
figurations defined by the variables m and n . The quantity m is the 
number of axial half-waves in the buckle pattern while n is the related 
number of full circumferential waves. For any single combination of stiff- 
ness values, an infinite number of load-wavelength combinations can be 
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possible* The critical load is the lowermost load which is just sufficient 
to hold the shell in a non-cylindrical deflected shape. Therefore the 
computation of a critical load must involve a mathematical minimisation. 

The nature of the problem is such that this minimization must be accomplished 
with respect to two independent wave-type parameters. The values m and n 
might be selected for this purpose although other choices, such as a and p, 
would be equally satisfactory. 

Since the analysis techniques of Volumes III [ 3 ] and IV t4] are 
based upon the minimization of equations (3-l) and (3-2), these methods 
only give classical small-deflection solutions. That is, the critical load 
predictions are obtained by locating a lowermost bifurcation point along 
the initially linear equilibrium path of a perfect orthotropic cylinder. 

Use of this approach raises some important questions as to the influences 
from initial imperfections and their interrelationship With the shape of 
the postbuckling equilibrium path. This matter is taken up in detail in 
Volume V Tl2] of this report. 
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SECTION 4 

OVERALL ANALYSIS PROCEDURE 

The methods presented in the several volumes of this report provide 
practical means for the analysis of instability in axially compressed 
circular cylinders having eccentric stringers and/or rings. In order to 
actually apply these methods, one should proceed as follows: 

(a) For cylinders having both stringers and rings, first compute 
the critical buckling stress for the curved isotropic skin 
panels which lie between the stiffeners. This value can be 
calculated by using the methods of Volume II [l3l. These 
methods directly provide design values that have been reduced 
from classical results in order to account for the effects of 
initial imperfections. Therefore, no further knock-down factors 
need be applied to these particular stress values. Although 
buckling of the isotropic skin panels is usually not catastrophic, 
the related critical stress values must be established in order 
to compute effective skin widths for use in the investigation 

of other possible modes. 

(b) Then determine the failing stresses for sections which lie 
between rings. For configurations which have no stringers, the 
monocoque cylinder curves given in the appendix of Volume IV [ 4 ] 
may be used for this purpose. These curves likewise provide 
design values which incorporate reductions from classical theory 
to account for the effects of initial imperfections. No further 
knock-down factors need be applied to these particular stress 
values. For configurations which include stringers, the methods 
of Volume III [ 3 ] should be used to obtain the failing stresses 
for sections lying between rings. It should be noted that these 
methods account for the possibility that crippling might occur. 
Note however that Volume III [ 3 ] only presents the results from 
classical theory and that, for design purposes, these values must 
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be reduced in accordance with the knock-down criteria of 
Volume V [12], Whether or not the configuration includes 
stringers, at this point in the investigation it is assumed 
that the rings experience no radial displacement. The validity 
of this assumption will be tested in the step which follows 
below. 

(c) By using the contents of Volume IV [4], the general instability- 
(see GLOSSARY) stress for the overall stiffened cylinder must 
then be determined. Here again, the methods of Volume IV [ 4 ] 
only give classical values and, for design purposes, these values 
must be reduced in accordance with the knock-down criteria of 
Volume V £ 12 ]. 

(d) For design purposes, it is then assumed that catastrophic collapse 
of the structure will occur at the lower of the two values obtained 
from (b) and (c) above. That is, depending upon the particular 
geometric proportions, the failure can occur either by the panel 
instability (see GLOSSARY) or the general instability (see 
GLOSSARY) mode. 
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